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Abstract. We consider the zero mass limit of a relativistic Thomas-Fermi- 
Weizsacker model of atoms and molecules. We find bounds for the critical 
nuclear charges that insure stability. 



1. Introduction 

The zero mass limit of the relativistic Thomas-Fermi- Weizsacker (henceforth 
ultrarelativistic TFW) energy functional for nuclei of charges zi > (which need 
not be integral) located at Ri, i = 1, . . . , K is defined by [4] 



(1) £(p) = a 2 J (Vp 1/3 ) 2 dx + b 2 / p 4/3 dx - I V{x)p{x) dx + D(p, p), 
Here p(x) > is the electron density, 



K 



(2) n*)=«E 

i=l '~ 

the electrostatic potential created by the nuclei 



(3) D(p,p) = ^J 



p(x)p(y) , , 

dx dy, 



-y\ 

the electronic repulsion, and a = e 2 /he « 1/137 is the fine structure constant. In 
units in which h = c = 1, the constants a 2 , and b 2 in (TT]) are given respectively by 

(4) a 2 = ^(3tt 2 ) 2 / 3 A 
and 

(5) b 2 = \{^ 2 f'\ 

In the non relativistic case, some emphasis has been placed on the question of an 
appropriate choice of the coefficient A of the gradient term connected to the kinetic 
energy. The non relativistic gradient correction A(Vp) 2 //3 was initially derived by 
Weizsacker [Tl] with a value A = 1, whereas the systematic gradient expansion by 
Kirznits [6] leads to A = 1/9. In the derivation of Tomishima and Yonei [13] , 
A = 1/5. Lieb [3 [9] showed, that the Weizsacker term introduces a z 2 correction 
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to the leading z 7 / 3 term of the non-relativistic ground state energy. Adapting the 
coefficient in front of the Weizsacker term such that the correction agrees with the 
leading z 2 correction of the non-relativistic quantum mechanics, the Scott correc- 
tion, lead Lieb to propose A = 0.185. In relativistic quantum mechanics, the leading 
energy correction remains unchanged whereas the Scott correction is smaller than 
in relativistic quantum mechanics [5], i.e., we cannot expect to have A the same 
value as in the TFW functional. However, we are not yet in a position to proceed 
with Lieb's strategy and to infer the coefficient of the gradient correction from this. 
In particular that would require showing that the massive equivalent of the function 
£ leaves indeed the leading energy contribution unchanged and the gradient term 
yields again a z 2 correction. 

Let us first consider the atomic case, i.e., the case K = 1, Z\ = z, Ri = 0. 
Because of simple scaling considerations, if we minimize the energy functional ([1]) 
over all functions p for which each of the terms in |T]) makes sense, the infimum 
of the energy functional is either zero or minus infinity. In the first case we say 
the atom is stable. Otherwise we say the atom is unstable. Our purpose here is to 
determine the range of values of the z^s for which the atom or molecule is stable. 
The following result holds in the atomic case (i.e., for K = 1, z\ = z, R\ = 0) pQ. 

Theorem 1.1. Let 

(6) £(» = a 2 J (Vp 1/3 ) 2 rfx + b 2 J pV*dx - J za-^dx + D(p, p), 

with D(p, p) given by f3|). Then 

-oo for z>%t + ^ 
for z < IF 

•> da 

where the infimum is taken over all nonnegative functions p(x), such that p S 
i 4 / 3 (R 3 ) ; Vp 1 / 3 S L 2 (R 3 ), and D{p,p) < oo. 



(7) 



inf £(p) = 



Remarks 

i) It follows from ([7]) that if z < 4a6/(3a) the atom is stable, whereas if z > 
4ab/(3a) + 77ra 3 /(66 3 ) the atom is unstable. The exact critical value of z (z c say) 
dividing the region of stability from the region of unstability is not known. However, 
it turns out that for the physical values of the constants the gap between the upper 
and lower bounds on z c is less than one, and therefore negligible (see the following 
remarks) . 

ii) For the physical values of a and b given by ^ and ([5]), the atom will be stable 
if z < vW^/a ~ 167.8a/A. Thus, if A = 1/9 (i.e., the value used by Kirznits, 
[TJ, [B] ) the atom is stable if z < 56. If A = 1/5 (i.e., the value used by Tomishima 
and Yonei, [13]) the atom is stable if z < 75. Finally, using the value of Lieb 

the atom is stable if z < 73. 

iii) As for the value of the gap, using the physical values of the constants, one gets 
77ra 3 /(36 6 ) = (7/12tt) v /3A 3 /2 < 0.021 < 1 for all the values of A considered above. 
Thus, the gap is negligible from the physical point of view. 

For the molecular case, i.e., when K > 1, and V is given by @ the following 
result was proven in [T]. 
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Theorem 1.2. Let 

(8) Si{p) = a 2 JiVp^fdx + b 2 J p 4 / 3 dx- J Vpdx + D(p,p), 
with V given by fl[) and D(p, p) given by |3J). Then 

(9) inf£(p)=0 tf£ = ££i*i<$£. 

where the infimum is taken over all nonnegative functions p{x), such that p G 
L 4 / 3 (M 3 ), Vp 1 / 3 e L 2 (R 3 ), and D{p,p) < oo. 

The above result is just a trivial extension of the atomic to the molecular case, 
and in some sense is the best possible when the interaction between the nuclei is 
not taken into account. In fact, if we neglect the nuclear interaction one can always 
think of the possibility of putting all the nuclear charges at the same point and 
reducing the molecular case to the atomic case, which explains the condition (J9]) 
on Z = J2i=i z i- The deficiencies of the above result are obvious. The goal of this 
paper is to have a result that yields stability for reasonable values of the nuclear 
charges. For that purpose the nucleus-nucleus interaction 

(10) U = a 



l<i<j<K 



\Ri — R 



plays a key role, because it prevents the possibility of putting the nuclear charges 
on top of each other. 

Our main result is the following theorem for the molecular case. 

Theorem 1.3. Let £(p) be given by 0j for functions p as in Theorem ] 1.21 with V 
given by Then, we have stability, i.e., 

inf£(p) + E/> 

provided 

4 a , 

(11) < z t < —bVT^x 

3a 

where x € (0, 1) is the root of 
(12) 



1-x b 4 /4^ 2 



3 J 27ra(4 + 9a 4 ) ' 



Remark For the physical values of a and b given by ((4]) and (O, and taken the 
physical value of the fine structure constant (i.e., a = 1/137) in (fT2|) . Theorem II. 31 
says that the molecule will be stable if each Zi < 55, if A = 1/9 (i.e., the value used 
by Kirznits, [BE]). If A = 1/5 (i.e., the value used by Tomishima and Yonei, [T3] ) 
the molecule is stable if each Zi < 74. Finally, using the value of Lieb [51 [5], the 
molecule is stable if each z% < 71. These bounds on each individual nuclear charge 
are almost the same as those embodied in the atomic case (i.e., the ones given in 
Theorem 11.11 above) . 

In the next section we give the proof of Theorem 11.31 
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2. Improved results on the stability of molecules 

In this section we prove Theorem II .31 Our proof relies in a modified uncertainty 
principle (see Theorem 12.11 below) which is of independent interest. As we men- 
tioned in the introduction, the nucleus-nucleus interaction plays a key role in the 
stability of molecules. As in |2J we may use the fact that the energy is separately 
concave in the nuclear charges and each charge Zi varies between and z. The 
minimum of a concave function is always on the boundary and hence the value od 
Zi wants to be either or z. If it is zero we have one nucleus less and if it is z then 
we are in the case we are considering. 

First we need some notation. We introduce the nearest neighbor, or Voronoi, 
cells [H] {Tj}f =1 defined by 

(13) T j = {x\\x-R j \<\x-R k \}. 

The boundary of Tj, dTj, consists of a finite number of planes. We also define 
the distance 

(14) Dj = dist (Rj , dTj ) = i min{ | R k - Rj | | j ? k} . 

One of the key ingredients we need in the sequel is an electrostatic inequality of 
Lieb and Yau [HJ [12] . Define a function <E> on R 3 with the aid of the Voronoi cells 
mentioned above. In the cell Tj, $ equals the electrostatic potential generated by 
all the nuclei except for the nucleus situated in Tj itself, i.e., for 16 Tj, 

K 

(15) = z^2 \x- Ril' 1 . 

i=l 

If v is any bounded Borel measure on R 3 (not necessarily positive) then 

(16) \f f \x-y\- 1 dv{x)du{y)- f $(x)dv(x) + U >\z 2 y^Dj\ 

We will also need a localization result for the kinetic energy which will allow us 
to control the Coulomb potential near each nuclei. This localization result for the 
UTFW model is given by Theorem 12.11 below . 

Theorem 2.1 (Modified uncertainty principle). For any smooth function f on the 
closed ball Br of radius R we have the estimate 

a 2 f \Vf(x)\ 2 dx + b 2 f f{xfdx>abf 

J Br J B ft J Bji 

Remark Notice that the factor 4/3 is best possible and it agrees with the sharp 
value given in Theorem 1.1. 

To prove this theorem, we need the following preliminary result. 

Lemma 2.1. Let u(r) be any smooth function with u(R) = 0. Then the following 
uncertainty principle holds 

\f (iudxD+uWxDIxDfixfdxiKiif \Wf(x)\ 2 dx)^ 2 (f u(\x\) 2 \x\ 2 f(x)Ux)^ 2 . 

J Br J Br J Br 
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f{xfdx. 
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There is equality if and only if 

f 



\J^[su(s)]ds + C 
for some constants C and A. 
Proof. Set 

gi(x) = u(\x\)xi 
where u is a smooth function with u{R) = and note that 



(3u(\x\)+u'(\x\)\x\)f(xfdx = Y, / f(x)[d j9j {x)]f{xfdx 

t j JB R 

= E/ f(x)d j (g j f)(x)dx-2j2 I f{xfg 3 {x)d 3 f{x)dx . 

j J Br j J Br 

Integrating the first term by parts yields 

W f(x)d j (g j f*)(x)dx = -J2 I d J f(x)g J (x)f(x) 2 dx+ [ f{x) 3 u(\x\)\x\dS(x), 

j JBr j JBr JdBn 

where the boundary term vanishes since u(R) = 0. Thus, 

f (3u(\x\)+u'(\x\)\x\)f(xfdx = -3j2 [ f{xfg j {x)d j f{x)dx. 

JB R j J B a 

Using Schwarz' inequality on the last term yields 

\ 1/2 / r \ 1/2 



<3 / \Vf(x)\ 2 dx\ / u(\x\) 2 \x\ 2 f(x) A dx 



(3u(\x\)+u'(\x\)\x\)f(x) 3 dx 

i 

Schwarz's inequality is an equality if and only if 

9jf = -\gj{x)f(x) 2 , 
which can easily be integrated and yields the stated function. □ 



Proof of Theorem \2.1\ To prove the theorem, pick 

U{r) =2 { r-R ] 
in the lemma which leads to the inequality 

3 



\x\ ~ 2R 



f{xYdx < 



\Vf(x)\'dx 



B r , 



1/2 



Bj 



{l-^) 2 f{xfdx 



1/2 



Next, applying the inequality between the arithmetic and geometric mean yields 

2 



[ \Vf(x)\ 2 dx + b 2 [ 

JBr J Br 



1 - 



R 



f(x) dx > ab 



B r 



4 2 
3~R ~ R 



f(x) 3 dx , 



from which the theorem follows. 



□ 
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In order to prove our main result, we consider the total energy, + U, and 
we split the J p 4 ! 3 term in two parts, i.e., {b\ + 6|) / p 4 ^ 3 - For later discussions, 
it is important to remark that the parameter 62 can be chosen arbitrarily in the 
interval (0,6). Then we use Theorem 12.11 with f 3 = p and Br — Bj, the largest 
ball inscribed in the corresponding Voronoi cell Tj, to get 

Z(P) + U> bjJ R3 p 4 / 3 dx-J R3 V(x)p(x) 

(17) + ab 2 Y,f=J B .(z^-%)p(x)dx + D(p,p) + U. 

In order to cancel the Coulomb singularity inside Bj we choose the parameter 

18 & 2 = T — • 

4 a 

The restrictions on 62 will give restrictions on z to insure stability. 
With the help of the Voronoi cells we now define, 

(19) W(x) ee + 

\x-Rj\ 

if x G Tj and \x — Rj\ > Dj, whereas, 

(20) W(x) ee *(x) + 

if x S Tj and \x — Rj\ < Dj. Now, if we restrict to values of z such that z < 
4a62/(3a), we can finally write 

(21) C(p) + tf >&(/>) + 6(p), 

with 

(22) ti{p)=bi f p 4/3 dx~a I W{x)p(x)+a [ <$>p{x)dx, 

JR3 JR3 Jm. 3 

and 
(23) 



= p)-a $(x)p(x) dx + U 

JM 3 

Using the Lieb-Yau electrostatic inequality (|16[) we have 
(24) fc(,)>a££;-L. 



3=1 J 

On the other hand, it is simple to estimate £i(p) from below, since it is simple 
to solve the variational principle inf p ^i(p). Thus we get, 

(25) ^^-r 4 {iSlj w -^ dx - 

From ([H]), (H5), and (|25J| we get 

> 3 ' 5 K 

(26) ^(p) + C />_^a ' 



Now, using PD|) we compute, 
(27) / (VF-$)|da;: 



2a6 2 \ 4 4 _ 3 64?ra 4 6^ 
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Since any Voronoi cell is contained in a half space, one calculates as in [KJ] that 

(28) / {W-<S>)\d x =( (^—) dx <^f. 
Finally, using the estimates (|2T|) and (|2"8"]l in (|2^|) , we get 

(29) ^P) + U >J M J2W.' 



3 = 1 J 

where 

(30) M = --a [3nz + 

If M > then £(p) + U > and the molecule is stable. Using the fact that 
b\ = b 2 — b\ together with (p~8|) the condition M > can be written solely in term 
of b\ as 

2 6 2 -6? /4\ 2 1 



(31) a z «— L - < - , , 

1 ; b\ ~ \3j 27r(4a + 9a 5 

Finally, in order to allow for the largest possible value of z (equivalently the largest 
possible value of 62), 61 must be chosen so that we have equality in (|3Tj) . Since the 
left side of (|3"Tj) is decreasing as a function of 61 in the allowed interval (0,6), we 
conclude the proof of our Theorem II. 31 



Remarks 

i) Note that we do not need a bound on the fine structure constant to ensure 
stability in this model. This is due to the absence of the exchange term. 

ii) If we set a = we find stability up to z = 71, when the parameter A = 0.185, 
up to z = 74 when A = 0.2 and z = 55 when A = 1/9. 
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